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Abstract 

The dimer problem arose in a thermodynamic study of diatomic molecules, 
and was abstracted into one of the most basic and natural problems in both statistical 
mechanics and combinatoric mathematics. Given a rectangular lattice of volume V in d 
dimensions, the dimer problem loosely speaking is to count the number of different ways 
dimers (dominoes) may be laid down in the lattice (without overlapping) to completely 
cover it. Each dimer covers two neighboring vertices. It is known that the number of 
such coverings is roughly e XdV for some constant A^ as V goes to infinity. Herein we 
present a mathematical argument for an asymptotic expansion for A^ in inverse powers 
of d, and the results of computer computations for the first few terms in the series. As a 
glaring challenge, we conjecture no one will compute the next term in the series, due to 
the requisite computer time and storage demands. 
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For simplicity and to be specific we consider a periodic cubical lattice A with the edge 
length even. Then the number of different coverings by dimers has as dominant behavior, 
in the V — > oo limit, 

e XdV (f) 

for a d-dimensional lattice, where the number of vertices we denote by both N and V. 
The asymptotic relation 

X d ~ \ ln(2d) - i , (2) 

was exactly presaged over 40 years ago by the results of [1], and this work was largely 
the basis for a complete proof of (2) a decade later in [2]. While in two dimensions A2 
is known exactly, [3], [4], in dimensions higher than two it is computationally difficult to 
get good bounds for A^. In a series of reports [5], [6], [7] we have found a mathematical 
argument for a full asymptotic expansion for A^ 

A,~i lB(M )-! + £l + ! + ... (3) 

and in [8] we presented the result of computer computations of the first three q leading 
to 

X d ln(2d) 1 + - + =r + ---. (4) 

2 v ; 2 8d 96 d 2 64 d 3 K J 

In the Appendix we show the consistency of (4) with known rigorous bounds on A^. We 
there also make some "intelligent guesses" for A 2 and A 3 based on the present theory, 
whose values further increase our confidence. We proceed to clarify and expand upon 
some of the developments in [5] - [8] . 

In d dimensions there are d 'kinds' of dimers, each kind oriented in one of the d 
lattice directions. Each kind of dimer may be 'located' at some place on the lattice. 
We generalize this situation as follows. A 'located tile' is a two element subset of the 
lattice. A 'tile' is an equivalence class of located tiles, with equivalence given by setting 
two subsets equivalent if one is a translation of the other. (A 'tile' is a generalization of 
a 'kind of dimer'.) Notice that tiles need not be connected. 

We consider tiles with a 'weighting', a function on tiles. We normalize the weightings 
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we consider, by requiring, if g is the weighting function, 

E 9(t) = 1/2 , (5) 

t 

the sum over all tiles, t. We let / be the weighting function given by 

if t is a dimer 

(6) 

otherwise 
that clearly satisfies the normalization condition, (5). 

A 'tiling' Tj of A is a set of two element subsets of A, 

Ti = {s\, s 2 , s l N / 2 } (?) 

where the s\ are disjoint. 

We now realize the sum over all possible dimer covers of A, the goal of our study, as 

(2d) N ' 2 Z (8) 

with 

z = e n ( 9 ) 

where the sum is over all tilings of A; the product over the f's selects those tilings in 
which all the tiles employed are dimers. The bar over a subset indicates the equivalence 
class of the subset as defined previously. 

We let fo be a constant function on tiles with value jf^j, to satisfy the normalization 
condition (5). We write 

f = fo + (f~fo) = fo + v, (10) 
/, f , and v all functions on tiles. Z becomes 

z = e n (/o +«(*«)) (ii) 

Z = Z + Zi + Z 2 + ... (12) 
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having expanded Z in powers of v. 



We note 

e NX d = z (13) 

Here A^ is understood a function of N, the usual A^ the infinity limit. Or equivalently 

\ d =l ln(2d) + ^lnZ (14) 

If one replaces Z by Z , the mean-field approximation in a natural nomenclature, one 
gets taking the infinity limit 

X d = \ ln(2d) -\ = \ ln(2d) + Jim 1 In Z (15) 

Z Z Z N^oo iV 

by an easy calculation. Thus equation (2) arises from the mean-field approximation, 
keeping only Z in (12). 

We return to (12) and introduce some convenient notations. 

Z = Z Z* (16) 

Z* = 1 + + + (17) 

Z* = Z t /Z (18) 

There is a natural factorization of Zi into a contribution from the factors of v in (11) 
which we call Z* and the factors of f in (11) which we call f3(N,i)Z so that 

Z; = (3(N,i)Z*. (19) 

A tedious calculation shows 

P(N,jN) ~ ^[(^Ka-^+i] (20) 
for large iV, a result we later need. 

We let Z* be Z* with /3(N,i) replaced by 1, 

Z* = 1 + Z{ + Z* 2 + ... (21) 
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where a detailed specification of Z* is given by 

a=l 

disjoint 

Now referring to [9] we may write a cluster expansion for Z* 

In Z* = ~7 Js (23) 



S! 

s 



J s = w(si)...w(s s )^(si,S2,--,s s ) (24) 



S1,S2,...,Sb 



where we may identify equation (2.5a) of [9] with Z*, and (23,), (24) with equation (2.7) 
of [9] . The located tiles in the sum of (24) are forced to overlap so they cannot be divided 
into two disjoint sets. ip' c is a numerical factor depending on the overlap pattern. To make 
our computations mathematically rigorous it will be necessary to study the convergence 
properties of sums such as in (23). At present this appears very difficult, and we by no 
means see yet a clear route to a proof. It is a challenging problem for the mathematical 
physicist. 

It is easy to show J\ — and it is proven in [6] that 

c c c 

with r > s/2. We also find it convenient to define 

NJ % = (1/i!) Ji (26) 



From (17), (19), (21), (23) one gets 

Z *= E P{N^ia l )J^---JZT — f (27) 

^ a x \ ■ ■ •a 9+1 \ 

We approximate the sum in (27) by its largest term, in the limit iV — > oo. If all the J's 
are positive this is a reasonable way to extract the dominant asymptotic limit. In [7] an 
argument is given that our results will hold even if some of the J's are negative. 
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We differentiate with respect to the ce's to find the largest term in (27), using (20) to 
treat the (3 factor and find the following equations, upon scaling c*j — > cti. 



d 

In a k = In Jfc + 



1 — 2 N Jo, , , 

ln(l — 2S loti) + S ioti 



(28) 



da k 
or 

a fc = J k e F * {a ' s) (29) 

where (28) and (29) define the F k . Extracting the term in the sum of (27) with the a k 
of (27) equal N times the a k of (29) we get 

n\ -S ch Fi + £ Ji e F i + 1-2 3 j in(l-2S i aA+Siai \ , N 

Z* ~ e I v y J . (30) 

Using (28) and (29) we can expand the exponent in (30) into a power series in the J's. 
From (14), (15), and (16) then we can find as a formal power series in the J's. Once 
the J's are known, one can extract a formal power series in the inverse powers of d (see 
(25)), and our formula (4). Everything rests on finding the J's ! 

All computer computations were done using maple in integer arithmetic. We note Jj 
is a function of dimension. We computed Ji, J2, Jq in dimensions 1,2, and 3 by com- 
puting the expression in (24). By a simple device one may replace the v's in (24) by /'s, 
and afterwards recover the original expression. Thus the most complicated computation 
of these involved placing down six dimers in three dimensions, in all possible ways where 
overlaps make it impossible to disconnect into two disjoint subsets of dimers. From the 
form of Ji as given in (25) we note this determines Ji,i = 1, ...,6 in all dimensions. To 
compute the c 4 term in (3) one would require J 7 and J 8 in arbitrary dimensions. This 
would entail computations in dimensions 1, 2, 3, and 4, the most complex setting down 
eight dimers in four dimensions. As another measure of difficulty of computing C4, we 
note it took approximately 10 seconds of computer time to compute c 2 , and two weeks 
to compute c 3 .... and we are using an exponential time algorithm. Good luck to the 
hearty soul attempting this computation. 

The computer results were first presented in [8]. The following table includes the 
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computations of the J's 





d=l 


d = 2 


d = 3 


Jl 











J2 


1/8 


1/16 


1/24 


J3 


1/12 


1/48 


1/108 


J 4 


- 3/64 


- 9/512 


- 5/576 


■h 


- 13/80 


- 23/1280 


- 11/2160 


J 6 


-19/192 


25/3072 


175/46656 



From these results there follows from (25) the following expressions for Jj 



Jl = 





(31) 


J 2 = 


1 1 

8 d 


(32) 


h = 


1 1 

12 rf2 


(33) 


J A = 


3 1 3 1 
_ 32 rf2 + 64 rf3 


(34) 


J 5 = 


11 3 1 
~ 8 d* ~ 80 # 


(35) 


■h = 


7 1 5 1 11 

48rf3 ~ 734^4 - q^5 


(36) 



Appendix 

We first study what rigorous upper and lower bounds on A^, for large d, have to say 
about our expansion (4). For an upper bound we have from eq. (9) and eq. (10) of [2] 
that 

x - < I ,n < M > -k+h Hil,d) + ^p- (37) 

Clearly an asymptotic expansion such as (3) with any value for the c's is consistent with 
(37), the 2 ln(d) dominating inverse powers of d asymptotically. 
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For a lower bound estimate one has from [10], and [11] Formula (5.4), [12] Equation 
(6.8) 

A d > - ((2d - l)ln(2d - 1) - (2d - 2)Zra(2d)) (38) 



from which follows 



1 1 °° 1 

A,>-fa(M)-- + g (M)( , :+1)(M) , (39) 

^M2d)-I + ii + li (40) 



For (3) to be consistent with (40) one must have 



ci > I (41) 



and if ci = 1/8 (as it does) then 



* > 5. (42) 

With ci = 1/8 and c 2 = 5/96 both inequalities, (41) and (42), are satisfied. A mini- 
mal test for the theory passed. The estimate (39) was pointed out to us in a personal 
communication by Shmuel Friedland. 

In extracting an asymptotic expansion in inverse powers of d from (28), (29), (30) 
we needed two J« values for each inverse power of d. (e.g. we needed Ji, J 6 to derive 
Ci, C3.) In [8] and [13] we introduced a formal parameter x (ultimately set equal to 1), 
replaced Jj by Jj x % ~ x in (28) - (30) and found an expansion in powers of x, this seemed 
a possible improvement. (The expansion in powers of x is formally somewhat like an 
expansion in powers of l/\fd.) 

We then viewed the expansion in powers of x, as an asymptotic expansion. Folio 
wing the notation of [8] and [13] we let Bk be the partial sum for A^ keeping powers of x 
through x k , then setting x to be one. We let B g and B g+ i be the two successive Bk with 
minimum value for \B g+1 — B g \. We set 

a = l( B 9 + B B+ i) (43) 
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and 

b = \B g -B g+1 \ (44) 

and took our estimate for A^ to be 

\ d = a±b (45) 

by the rule of thumb I learned as an undergraduate physics major. For d = 2 and d = 3 
this yielded 

A 2 = .296 ± .007 (46) 



and 

A2 is known exactly as 
by [3] and [4], and 
by [15], an improvement on [14]. 



A 3 = .453 ± .001. (47) 
A 2 = .2915.... (48) 
.440075 < A 3 < .457547 (49) 



Viewing the consistency of (46) with (48), of (47) with (49), and the validity of 
(41), and (42); we may associate a casual probability of ^, ^, |, \ for each of these four 
consistencies to occur at random. So we might say the odds are 6000 to one our theory 
is correct. While we may have presented this argument in a frivolous manner, some 
encouraging points are being made. (Using the consistency of (46) - (49) as an argument 
for the theory was suggested to me by Mihai Ciucu.) 

Acknowledgment. I would like to thank Shmuel Friedland, Gordon Slade and Nathan 
Clisby for support along the way. 
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